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A NOTE ON SECONDAEY FLOW IN ROTATING RADIAL CHANNELS^ 

By Jambs J. Kbambb and John D. Stanitz 


SUMMARY 

A general vector differential eguation for the vorticity com- 
ponent parallel to a streamline is derived for steady, nonviscous, 
and incompressible flow in a rotating system. This equation 
is then simplified by restricting it to rotating radial channels 
and by making further simplifying assumptions. The sim- 
plified equation is used to solve for the secondary vorticity, the 
vorticity component parallel to the streamline, in three special 
cases involving different streamtube geometries; the results are 
presented in a series of figures. The secondary vorticity is 
shown to decrease with decreased absolute angular velocity of 
the fluid, decreased inlet total-pressure gradient, decreased 
length of relative flow path, and increased relative velocity. 

INTRODUCTION 

When a fluid with a total-pressure gradient normal to the 
osculating plane of the streamline follows a curved path, 
velocities are induced in the surfaces normal to the irrota- 
tional or “through” flow direction. These induced velocities 
or secondary flows are recognized as the source of several 
types of loss in turbomachinery. For example, losses may 
occur because of (1) boundary-layer separation caused by 
the transfer of low-energy air to regions of decelerating flow, 
(2) improper angles of attack in the blade rows of compres- 
sors and turbines, and (3) viscous dissipation and unrecover- 
ability of the kinetic enei^ involved in the secondary flow. 

Because secondary flows are an important source of loss, 
considerable analytic work has been done on this phenomenon 
in stationary curved channels by Squire and Winter (ref. 1), 
Hawthorne (ref. 2), Exonauer (ref. 3), and others. The 
general method used in these analyses is to compute the 
vorticity component parallel to the actual streamline, the 
secondary vorticity, caused by a variation in hflet total 
pressure across the passage normal to the osculating plane 
of the streamline. This secondary vorticity is computed 
because it is an indication of the magnitude of the secondary 
flow. In addition, Exonauer (ref. 3), using a similar method, 
has computed this vorticity component in an axial-flow 
rotating channel such as that in an axial-flow compressor. 
However, no work has been done on radial-flow rotating 
channels, such as in a centrifugal-type impeller. 

This report extends the work of Hawthorne and the others 
to a rotating radial channel. The results of the analysis are 


indications' of the qualitative trends of variables which 
affect the secondary vorticity and therefore the secondary 
flow. The purpose thus is not to obtain exact data but 
rather to obtain some insight into the problem of secondary 
flows in rotating radial channels. It is hoped that deter- 
mining the relative importance of some of the variables 
involved will lead to a better understanding of the flow in 
centrifugal-type impellers and similar flow machinery. 
Calculations are carried out to show the effects of these 
variables in several streamtubes with various geometries. 

METHOD OF ANALYSIS 

In this analysis the distribution along a streamline of the 
secondary vorticity (component of the vorticity parallel to 
the stre amlin e) is computed for rotating radial channels. 

PEKLIMINABT CONSmEBATIONS 

For steady three-dimensional flow of a nonviscous fluid 
with nonuniform total-pressure distribution through a 
stationary channel, the streamlines and vortex lines lie on 
surfaces of constant total pressure (ref. 4, p. 244) called 
Bernoulli surfaces. The total-pressme gradient VP is neces- 
sarily normal to the Bernoulli surfaces, and in general these 
. surfaces are curved in space. The vorticity vector f, which 
is tangent to the vortex line, is defined by 

F=vxF 

where Vis the absolute velocity vector tangent to the stream- 
line. From the equation of motion the vectors E and 
VP are related by 

fx7=vp 

Since VP is normal to the Bernoulli sinface, this equation 
shows that the streamlines and vortex lines lie on Bernoulli 
surfaces. In general, the vorticity vector f has components 
tangent and normal to the absolute velocity vector V, f cos 
and f sin 1 ^-, respectively, where i/ is the angle between the 
vortex line and the stroamline (see fig. 1). The vorticity 
component tangent to is the secondary vorticity, and it is 
this component of vorticity that is an indication of the 
magnitude of secondary flow. Between any two points on 
the Bernoulli surface there exists a geodesic or path of 
minimum distance. In general, the streamline is not a path 
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of miniimim distance on the Bernoulli surface, that is, t^ 
streamline is not a geodesic; therefore the principal normal N 
of the streamline is inclined at an angle <p (fig. 1) from the 
normal to the Bernoulli surface VP. 

If the streamline is not a geodesic on the Bernoulli surface, 
the secondary vorticity increases in the direction of the flow 
(ref. 2). This condition exists if is not zero, that is, if VP 
has a component normal to the osculating plane of the 
streamline, the plane _containing the velocity vector V and 
the principal normal N. If the streamline is a geodesic, VP 
lies in the osculating plane and the secondary vorticity does 
not change along the streamline. 



The general assumptions and method used in the develoj)- 
ment of the equations for secondary vorticity in a rotating 
radial channel are outlined in the following sections. 

General assumptions. — ^The flow is assumed to be steady, 
incompressible, and nonviscous. Although the fluid is 
assumed to be nonviscous in the region of the solution, a 
gradient in inlet total pressure normal to the flow planes is 
assiuned that woifld have been caused by the generation of a 
thick boundary layer by viscous forces upstream of the 
region of solution, that is, upstream of the rotating channel. 
If it is assumed, as will be done herein, that the boundary- 
layer thickness is large compared with the passage height 
6. e., viscous shear forces are low), then the fluid can be 
assumed nonviscous and the vorticity calculated from the 
equation of motion for an ideal fluid. This assumption is 
also made in references 1 to 3. 

The flow entering the passage is assumed to be parallel 
and uniform except for the thick boundary layer which 
results in a total-pressure gradient normal to the osculating 
plane of the relative streamline. In a centrifugal-type 
impeller, this would be the boundary layer on the inner and 
outer shrouds. Thus, initially, the vorticity component 
parallel to the streamline is zero. A component of vorticity 


parallel to the streamline is generated when the streamline 
turns from its initial direction or is rotated m a rotating 
system. 

Outline of method. — In order to evaluate the secondary 
vorticity along a streamline, the vorticity component 
parallel to a streamline is expressed as a fimction of the 
relative total vorticity and the relative velocity. This 
expression is simplified using the fact that the divergence of 
the total vorticity is zero. The total relative vorticity is 
eliminated from the equation by means of the equation of 
motion. The equation of continuity is used to express the 
gradient along a streamline of the magnitude of the vorticity 
component parallel to a streamline divided by the relative 
velocity as a function of the relative velocity, the angular 
velocity of the rotating system, and the gradient of the 
Bemoulh constant in the rotating system. This is the 
general vector differential equation for secondary vorticity 
in a rotating system. This general vector differential 
equation is then simplified by restricting it to rotating radial 
channels and by making further simplifying assumptions. 
In order to solve for the secondary vorticity by means of 
this equation, it is necessary to know the streamline pattern. 
As the actual streamline configuration is, of course, unkno^vn 
at the outset of the problem, it is necessary to assume a 
realistic streamline pattern which from general considerations 
would probably not differ greatly from the actual pattern in 
order to obtain meaningful results. This assumed stream- 
line configuration is based on a potential flow solution or 
some other approximation, and a first approximation to the 
secondary vorticity is obtained by means of the simplified 
differential equation for secondary vorticity. The secondary 
vorticity is solved for in several rotating radial streamtubes 
with various geometries using the simplified differential 
equation. 

GENEBAn VECTOR DUTEBENTIAL EQUATION FOR SECONDARY 
VORTICITY 

The various steps outlined in the development of the gen- 
eral vector differential equation for secondary vorticity are 
now presented in detail. 

Vorticity vector. — The total relative vorticity is finst ex- 
pressed as the sum of its components parallel and normal to 
the streamline. Let W_be the velocity vector relative to 
the rotatog system and fi, the relative vorticity vector equal 
to VXW. (All symbols are defined in appendix A and are 
expressed in terms of dimensionless ratios.) The dimension- 
less relative velocity W is a ratio of the dimensional relative 
velocity to some characteristic speed of the channel ur,. The 
angular velocity of the system is co and r, is a reference dis- 
tance from the^axis of rotation to some point in the channel. 
The vorticity n and the magnitude of its component { parallel 
to a streamline are dimensionless, having been divided by to. 

It can_^e seen from the expansion of the triple vector 
product (WXQ)X^that 

w-w w-w 
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From the definitions of the dot and cross products it follows 
that the terms on the right are equal to the components of 
vorticity parallel and normal, respectively, to a streamline. 
Since { is the magnitude of the component of the vorticity 
parallel to a streamline, the preceding equation can be 
written 


^7 1 C^Xo)XPT 
+ WW 


( 1 ) 


where the absence of the bar indicates a scalar quantity. 
Because the divergence of the curl of any vector is zero, ft 
can be eliminated from the left side of equation (1) as follows: 


V-ft = V-(vxW) =0 
so that equation (1) becomes 


V-^W=V 


T7X(WX^ 
WW 


( 2 ) 


In order to simplify equation (2) the equations of motion 
and continuity are introduced. 

Equation of motion. — The equation of motion can be 
written (ref. 5, eq. (15)) 

WXft = Vff+^^ (3) 

CO 

where H is the Bernoulli constant relative to a rotating sys- 
tem, defined symbolically as 



Figubb 2. — Vector diagram for rotating ohanneL 


and, from equation (5) , 

W-W 


General differential equation for secondary vortioity. — 
Equations, (6) and (7) are now substituted into equation (2), 
and the following differential equation results: 


H= 


. vIp , TP 

■(wf,)'"^ 2 2 


(4) 


ir-v^v.^5+F.v 


-2S-F 


( 8 ) 


In equation (4), p/p is the ratio of the static pressure to the 
fluid mass density, made dimensionless by dividing by 
(orO*, and B is the dimensionless distance in the radial 
direction expressed as a ratio of r„ as are aU linear dimen- 
sions. The quantity H is analogous to the total pressure P 
in the stationary channel in that the relative streamlines lie 
on surfaces of constant H in the rotating case (fig. 2) just 
as the streamlines lay on surfaces of constant P in the sta- 
tionary case (fig. 1) . From equation (3) it can be seen that 

WX (WXft) = W X 

0 ) 

=FXVfl'+2lP--2 — F (5) 

CO 0) 

By means of equation (6), the right side of equation (2) can 
be expressed as a function of W, w/w, and VH. 

Equation, of continuity. — In order to solve for the com- 
ponent of the gradient of ^/W parallel to the streamline, the 
equation of continuity is introduced. This equation states 

V-F=0 

Hence, 

( 6 ) 


The first term on the right of equation (8) represents the 
secondary vorticity growth caused by two effects: (1) the 
turning of the streamline caused by the rotation of the 
system, and (2) the turning of the streamline caused by its 
curvature in the relative system. In order to clarify the 
physical significance of equation (8), the first term on the 
right is expanded as shown in appendix B to yield 


WXVH 


2 


CO 


2k,|VH1 

w 


(9) 


where is the geodesic curvature of the streamline on the 
Bernoulli surface. Combining equations (8) and (9) results 
in the following general vector differential equation for 
secondary vorticity: 


Tins vector differential equation relates the secondary 
vorticity { at a point along the streamtube to the angular 
velocity of the system w, the relative velocity of the fluid W, 
the geodesic curvature of the streamtube Xg, and the gradi- 
ent of the Bernoulli constant VH. 
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Physical interpreta,tioii of general vector differential 
equation. — ^The first term on the right of equation (10) is 
influenced by the component of VH parallel to the axis of 
rotation, and the second term is influenced by the component 
of VH normal to the principal normal of the streamline. 
The third term on the right of equation (10) represents the 
variation of the vorticity component parallel to the stream- 
line caused by the rotation of the system, irrespective of any 
variation of H. This type of secondary motion" is discussed 
in reference 6 for an asdal-flow-type channel. 

From the first two terms on the right of equation (10), it 
can be seen that in a rotating radial channel such as in a 
centrifugal impeller, in which the osculating planes of the 
relative streamlines are approximately normal to the axis 
of rotation, the only component of the gradient of H which 
affects the generation of secondary vorticity is that parallel 
to the axis of rotation. This corresponds to the component 
of the gradient of H from the inner to the outer shroud. Thus, 
the boimdary layer on the blade surfaces with the concomi- 
tant gradient of H normal to the blade surfaces does not affect 
the generation of secondary vorticity. 

The vorticity component parallel to the streamline induced 
by the rotation of the system is equal to — 2w cos a, where a 
is the angle between the streamline and the direction of the 
axis of rotation (see fig. 2). The applications in this report 
are to rotating radial channels so that this term is zero 
because a is 90°. If the channel is stationary, that is, w is 
equal to zero, equation (10) becomes 

W-V 

where it is understood that the variables are in dimensional 
form because w is equal to zero. This is the same equation 
as that obtained by Hawthorne in reference 2 for a stationary 
curved channel. 

APPUCATION OF GENEEAL VECTOK DtFFEHENTIAI. EQUATION TO 
HOTATING BADIAL CHANNELS 

In this section the general vector differential equation for 
secondary vorticity developed in the preceding section is 
applied to rotating radial channels similar to those in 
centrifugal-type impellers, an example of which is shown in 
figure 3, and further simplifying assumptions are made. 
Numerical examples are solved involving streamtubes "with 
various geometries which follow logarithmic spiral paths, 
including a straight radial path as a special case. A 
logarithmic spiral path is similar to the streamlines in a 
conventional centrifugal impeller and an analytic expression 
is known for such a path. 

Further assumptions. — ^A right-handed orthogonal curvi- 
linear coordinate system with coordinates «i, %i, and Uz is 
chosen. The unit vectors in the Ui-, Ur, and ■Wj-directions 
are 1, ~j, and respectively. The reference point denoted 
by the subscript t m the general analysis shall refer to the 
tip of the rotating radial channel so that aU linear dimensions 
are expressed as ratios of the tip radius and all velocities are 
expressed as ratios of the tip speed wr,. 


% _2Kt\VE\ 

W~ W 


Bernoulli surfaces 
approximately 
• normal to oxis 




Fiqubb 3. Front and side views of centrifugal-type impeller. 



Figueb 4. — Flow planes for simplified flow in rotating radial channel. 


The flow surfaces on which the relative streamlines lie are 
Bernoulli surfaces. These sinfaces are parallel to the UiUr 
surfaces and are assumed to be planes normal to the axis of 
rotation (fig. 4). As a result of this assumption, 


and 


o}=wk 


( 11 ) 




( 12 ) 


where the C/, for n=l, 2, 3 are the square roots of the 
components of the metric tensor associated with the orthog- 
onal curvilinear coordinate system employed and are given 

by 
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( 13 ) 


in wMch V is expressed in the coordinates in which the 
are given. From this assumption it follows that the 
Bernoulli smfaces are planes coincident with the osculating 
plane of the relative streamline throughout the channel and 


tp is +90° or —90° when 


(1.'^ 
\jJz bUiy 


is positive or negative, 


respectively. This assumption, which is usual in theoretical 
studies of secondary flow, imphes that the secondary flows 
are small perturbations on the primary flow. Hence, if tiie 
secondary flows are large, the results of this analysis are 
qualitative rather than quantitative. 

The orthogonal cmwilinear coordinates in the ■ai'Uj-plane 
are selected so that lines of constant Ut are streamlines 
Hence, 

( 14 ) 


It is not necessary for the application of the method that the 
coordinate system be so chosen but is merely a matter of 
convenience. As stated previously, it is necessary to assume 
a streamline pattern in order to solve for the secondary 
vorticity. The accuracy of the solution for secondary 
vorticity is therefore dependent upon the accuracy of the 
streamline approximation. This approximation might be a 
two-dimensional potential flow solution, a mean flow path 
through the channel, potential flow through the stationary 
channel, and so forth. 

Although the Bernoulli surfaces are assumed to be essen- 
tially parallel planes, the gradient of the Bernoulli constant 
is oUowed to vary inversely as the assumed streamtube 
thickness ratio r in order to give a closer approximation to 
conditions in an actual passage. The streamtube thickness 
ratio T is the ratio of the streamtube thickness in the 
'Us-direction at any point to its value at the inlet of the 
rotating radial channel. Thus, 


j_ ag_i /J^ dg \ 

Uz ^Uz T ^^ 3 / i 


( 16 ) 


Simplified differential equation for secondary vorticity 
in a rotating radial channel. — ^Under the assumptions and 
definitions given by equations (11), (12), (14), and (15), the 
differential equation (10) for secondary vorticity in a rotating 
radial channel becomes 


vector. It can be seen from equation (16) that the differen- 
tial change in secondary circiflation parameter d{^/W) is 


which is the Bemmflli 


directly proportional to - ^ 

constant gradient. Also, the differential change in the 
secondary circulation parameter is proportional to the sum 
(l+TFic) and to UiduilW. In the sum (1 +TFk), the term Wk 
is the angular velocity about the axis of rotation of the fluid 
relative to the rotating system, and the term 1 is the angular 
velocity of the system itself. Thus, the sum is equal to the 
absolute angular velocity of the fluid particle. Because the 
factor UiduilW is the differential element of time required 
for the fluid particle to pass through the differential channel 


length Uidui, the product (1 +TFk) 


Uidtii 

~w~ 


is the differential 


absolute angle through which the fluid particle is turned. 
I n addition, d(^IW) varies inversely with TP so that large 
values of di^jW) will occur for small values of W. 

If the velocity W is constant, the secondary vorticity f 
varies linearly with the secondary circulation parameter f/T7. 
If the flow is accelerating or decelerating for a given differen- 
tial change in f/TF, the differential change in |fl is greater or 
less, respectively, than that for constant W, as can be seen 
from the second term on the right in the foUowing expansion 
of di^jW): 

d^=.W d dW (17) 


However, becaxise the differential change in secondary cir- 
culation parameter di^jW) varies inversely with W (see eq. 

(16) ), the tendency of the first term on the right of equation 

(17) is to produce an effect opposite to that of the second 
term; that is, accelerating and decelerating flows tend to 
decrease and increase, respectively, |f|. Thus, there is a 
double effect involved so that it is not possible, in general, to 
predict the exact effect of accelerating or decelerating flow. 

logarithmic spiral coordinates. — A. logarithmic spiral 
curvilinear coordinate system is specified in the relative 
flow planes and a rectilinear coordinate system normal to 
the flow planes. Hence, the equations of the coordinates are 

■Ui=ln tan 

-if2=0— tan /3 In jB > (18) 

Uz=Z ) 



or along a streamline. 




Uidui 

W 


( 16 ) 


The term f/TF is proportional to the secondary circulation, 
that is, the circulation associated with the secondary 
vorticity, because f is the secondary vorticity and 1/W is 
proportional to the streamtube cross-sectional area normal 
to the flow and therefore normal to the secondary vorticity 


where B, 6, and Z are the usual cylindrical coordinates 
fo rming a right-handed system and |3 is the streamtube angle, 
that is, the angle between the radial and •Wi-directions, posi- 
tive when in the direction indicated by positive a according 
to the right-hand rule (see fig. 4). For logarithmic spiral 
streamtubes, is constant along the streamtube. 

The curvature x of a logarithmic spiral streamtube is 
sin P/R, and the differential streamtube length is sec /3 dB. 
Hence, along a logarithmic spiral streamline equation (16) 
becomes 
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The following equation results from integrating equation (19) 
along a relative streamline and noting that is equal to 
zero because the flow was assmned to be initially parallel 
and uniform except for a total-pressure gradient normal to 
the osciflating plane of the streamline: 

^=2Wr sec /3 ^ (h ^ (20) 

where ® is an arbitrary point along the streamline. This 
equation can be used to compute the secondary vorticity at 
any point z along the streamline from the known geometric 
and flow conditions. 

Case I. Constant t and WR . — ^In Case I, 

( 21 ) 

and the streamtube geometry is similar to the channel in a 
centrifugal impeller with constant blade height and blade 
spacing. In this case the central angle subtended by the 
streamtube width Ad is constant. Because t is constant, 
bHfdZ is constant along the streamline. Integrating 
equation (20) using W as determined by equation (21) 
yields 

Equation (22) gives the secondary vorticity $ at any radius 
B for Case I. It can be seen that $ is directly proportional 
to (bfl/bZ)t and that the change in | between any two radii 
is dependent on the value of i2«. 

Several special cases can be considered with this particular 
georhetry. If is equal to zero, equation (22) becomes 


and the streamtube geometry is s imil ar to the channel in a 
centrifugal impeller with constant blade spacing {A6 is con- 
stant) and constant flow area (blade height varies inversely 
with radius). From equations (15) and (26), bHjbZ is 
directly proportional to the radius R. Integrating equation 
(20) using equation (25) for t yields 

WB^t \w (^*— ^ (12— i2<)J (26) 


Equation (26) is the expression for secondary vorticity at 
any point in the channel for Case II. The change in { 
between any two radii is inversely proportional to the value 
of Rj. 

From equation (26) the secondary vorticity is zero at any 
point in the channel when 


14 


W sin /3 



0 


( 27 ) 


It is interesting to note from equation (27) that the second- 
ary vorticity is zero at a point when the absolute angular 
velocity (left side of eq. (27)) at the mean radius between the 
inlet and that point is zero. In the previous discussion of 
equation (16), it was noted that if IF is constant, the differ- 
ential change in secondary vorticity at a point is zero if the 
absolute angular velocity of the fluid was zero at that point. 

Case m. Constant W and t. — I n Case IH both the 
relative velocity W and the streamtube height ratio r are 
constants so that the streamtube geometry is similar to the 
channel in a centrifugal impeller with constant blade height 
and constant flow area (blade spacing A9 varies inversely 
with radius). Integrating equation (20), having noted that 
W and T are constants, yields 





\-R sin 




( 23 ) 




If the streamtube walls are straight radial lines, /3 is zero and 
equation (23) becomes 



Equation (28) expresses the variation in { along a streamline 
for Case HE. It can be seen that the change in { between 
any two radii is independent of the value of The condi- 
tion for zero ? at a point is 


Also, if the channel is stationary, <o is zero and equation (23) 
becomes 

„ l/'bH\ _tan /3 
^iK-dzJr W 

AH the quantities in the above equation must bo expressed 
in their dimensional form because to is equal to zero. 

The condition for ? to be zero at some radius R is given by 

^^^4-sin p=0 (24) 


as indicated by equation (22) . 

Case n. Constant tB and W . — In Case II, 

^~R 


(.25) 


, , TF sin /3 

^^(B-Bi) 




(29) 


The average absolute angular velocity between the inlet and 
radius R is given by 


B 


1 I 11^ sin ^ , IF sin d 1 B 


so that equation (28) indicates that the secondary vorticity 
is zero at a point when the average absolute angular velocity 
(left side of eq. (29)) is zero. This case corresponds closely 
in geometry to that considered by Hawthorne and others in 
stationary elbows. The condition for zero f in Case IH, 
namely, that the average absolute angular velocity of the 
fluid is zero, is analogous to the condition of zero turning 
angle for zero secondary vorticity in a stationary elbow. 
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Summary of equations. — summary of the geometric 
and flow conditions as well as the final equations for second- 
ary vorticity and the condition for zero f in each case is 
presented in the foUowing table: 



Variable 

Case I 

Case II 

Case m • 

Specified 

conditions 

A0 

Constant 

Constant 

(A9)t 

R 

T 

Constant 

Ri 

R 

Constant 

i)H 

Constant 

R/bH\ 

RA^zJ, 

Constant 

Consequences 
in flow 

W 

Wt 

R 

Constant 

Constant 


Eq. (22) 

Eq. (20) 

Eq. (28) 

Condition 
for {= 0 

Eq. (24) 

Eq. (27) 

Eq. (29) 


RESULTS AND DISCUSSION 

Standard conditions. — The final equations express the 
secondary vorticity { as a fimction of (bHfbZ){, Wt, -K, 
and /3. As has already been noted, the secondary vorticity 
is directly proportional to (bHfdZ)i. The effects of the 
other variables are presented in a series of figures in which 
the variation of f/(djS/dZ)< with two variables is shown with 
the other variables given their values at standard conditions. 
These standard conditions are: 

Wt 0.4 

/3 0 

Ei 0.4 

R 1.0 

Case I. — The variation of the secondary vorticity para- 
meter with radius ratio R for various values of 

the inlet radius Rt and the streamtube angle )3 is shown in 
figure 6 for Case I. In figure 5 (a) the effect of the inlet 
radius on the growth of the secondary vorticity in passing 
through the channel is shown. The variation of the second- 

ary vorticity parameter with inlet radius is small for this 
case, so the choice of the value of the inlet radius for the 
standard conditions will not affect the presentation of the 
results to any considerable extent. 

From the curve of f/(c)fl/5Z)< against R for various values 
of streamtube angle /3 in figure 6 (b), it can be noticed that 
the change in secondary vorticity at the tip radius for a 
given change in streamtube angle decreases as the streamtube 
angle decreases. This phenomenon is caused by the fact 
that although the absolute angiilar velocity component 
decreases as /3 decreases from 40° to —40°, the length of path 
through the channel at first decreases as /3 changes from 40° 
to 0 but then increases as /3 varies from 0 to —40°. As 
discussed following equation (16), the change in secondary 
vorticity varies directly with both the absolute angiilar 
velocity and the path length so that when /3 decreases from 0 


the effects partially cancel each other. For j3 equal to —40°, 
the secondary vorticity parameter becomes negative over a 
range of R following Fj because the absolute angular velocity 
is initially negative. 



(b) Effect of p. 

Figveb 6. — Variation of H(bnfd2^t with R for Case I. Standard 
conditions unless otherwise noted. 

Figure 6 (a) presents the variation in the secondary vor- 
ticity parameter at the tip with tip velocity 

ratio Wt for several values of /3. The effect of Wt is pro- 
nounced, as would be expected from equation (16). For W» 
equal to 0.2 and equal to —40°, is greater than for Wt 
equal to 0.2 and equal to —20° or 0, because the elongation 
of the flow path more than balances the effect of reduced 
absolute angular velocity of the fluid. 

Cases n and in. — The variation of f,/(dH/2>Z)j with Wt 
for several values of /3 is shown in figure 6 (b) for Case II. 
As in Case I, the values of ?t/(dS'/dZ)i for Wt equal to 0.2 
and equal to —40° are greater than the corresponding 
values for equal to 0 and —20°. The same reason men- 
tioned under Case I apphes here. Secondary vorticities 
are, in general, much greater in this case than in Case I 
because for the seme value of the gradient of the Bemoulh 
constant bHfdZ at the inlet, the value of ^HfdZ at any 
other point in the channel is greater in Case II than in Case 
I. Again the marked increase in secondary vorticity with 
decreasing Wt can be noted. Figure 6 (c) presents the same 
information for Case m as is shown in figures 6 (a) and 
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(a) Case I. 

Figubb 6. — Variation of {(/(di3/dZ)< with TFj. Rt equal to 0.4. 

(b) for Cases I and ]1, respectively. Again the secondary 
vorticity parameter increases rapidly for small values of Wt. 

This same trend is evidenced in the experimental results 
reported by Spannhake in reference 7 (pp. 145-159). The 
apparatus used for the experiment was an S-shaped tube 
that was rotated about an axis as shown in figure 7. The 
fluid, water, enters and leaves in the axial direction along 
the axis of rotation. A loss coefficient X based on the total- 
pressure drop from inlet to exit was computed. A range of 
ratios of circumferential velocity of the pipe to through 
flow velocity of the water, corresponding to 1/11^, from 0 to 
30 was used. The loss coefficient for flow through this 
apparatus shows the same trend, that is, rapid increase as 
1/W increases, as that of the secondary vorticity in the 
theoretical case. Since the experimental work is only some- 
what s imil ar to the theoretical work, for in a channel such 
as the rotating pipe the streamline configuration is much 


more complex, and since viscous dissipation, not considered 
in the theoretical work, is an important factor, the only 
conclusion which can be drawn is that the trends are similar. 

Comparison of Cases I, II, and IH. — A. comparison of the 
variation of with W, is shown by the solid lines 

in figure 8 for the three cases. It is immediately apparent 
that the secondary vorticity for Case I is less for all values 
of Wt than it is for Cases II and HI. If the secondary 
vorticity parameter {,/(dSydZ)j were plotted against inlet 
velocity Wt, the curves for Cases II and III would be the 
same since the velocity is constant throughout the stream- 
tube. The curve for Case I is shown as the dashed line in 
figure 8 for which the abscissa is Wt. It is seen that the 
curve for Case I lies above the curves for Cases II and III. 
As discussed following equation (16), the effect of decelerat- 
ing flow is to increase the value of d(^/W), because of the 
consequent smaller values of W, and at the same time to 



(b) Case II. 

Figubb 6. — Continued. Variation of with TV,. Bt equal 

to 0.4. 
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(o) Case in. 

FiauRB 6. — Concluded. ' Variation of with TV,. R{ equal 

to 0.4. 

decrease for a given value of d{HW). Thus it appears 
that of these two effects, that of increasing d(J^jW), and con- 
sequently df, predominates because ? is greater for Case I 
than for Cases H and lU when all three cases have the 
same inlet velocity, that is, when the velocity is less for 
Case I than for Cases II and 111 at all points in the stream- 
tube other than the inlet. When all three cases have the 
same velocity at the tip, as shown by the solid lines in figure 
8, the velocity is greater for Case I than for Cases II and 
m at aU points in the channel other than the tip. Thus 
is smaller and also the flow is decelerating so that, 
as can be seen from equation (17), these two effects com- 
pound, with the result that the secondary vorticity param- 
eter is smaller for Case I than for Cases II and DI. 




Figtjbb 7. — Variation^ of loss coefficient X with velocity ratio 1/TF 
measured in rotating S-ahaped tube as reported in reference 7. 
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The secondary vorticities for Case II are larger than those 
for Cases I and III because if (dS/bZ)* is the same for aU. 
three cases, bHfdZ at any other point is greater for Case II 
than for Cases I and IH. The secondary vorticity, being 
directly proportional to bH/bZ, is therefore greater for Case 
n. All three curves show the same pronounced increase in 
secondary vorticity as TF, decreases. 

SUMMARY OF RESULTS AND CONCLUSIONS 

A ‘general vector differential equation for the vorticity 
component parallel to a streamline in a rotating system is 
derived. This equation indicates that at a given point along 
the streamtube the secondary vorticity f is a function of the 
angular velocity of the system S, the relative velocity of the 
fluid W, the geodesic curvature of the streamtube k„ and the 
gradient of the Bernoulli constant VH. This equation is 
then simplified by restricting it to rotating radial channels 
and by making further simplifying assumptions. This sim- 
plified equation was used to solve for the secondary vorticity 
in three special cases involving different streamtube geome- 
tries, and the results are presented in a series of figures. 
These figures and the equations from which they were 
obtained indicate that: 

1. The differential change in secondary circulation param- 
eter f/lF at a point in a rotating radial channel is directly 
proportional to the absolute angular velocity of the fluid, 
the gradient bHfbZ of the Bernoulli constant, and the differ- 
ential element of time required for the fluid particle to pass 


through the differential channel length. The differential 
change in |/1F also varies inversely with the square of the 
relative velocity. 

2. The secondary vorticities are lowest for Case I ■with 
constant streamtube thickness ratio t and constant central 
angle subtended by streamtube width Ad because of the higher 
fluid velocities upstream of the channel tip. 

3. The secondary vorticities are highest for Case II with 
constant tB (where B is radius ratio) and W because the 
decreasing streamtube thickness downstream of the channel 
inlet results in increased bHfbZ. 

4. The secondary vorticity increased rapidly in all cases 
as the relative velocity at the tip decreased because of the 
inverse variation of df with the relative velocity. 

5. The decrease in absolute angular velocity of the fluid 
caused by a decrease in the streamtube oqgle /3 caused the 
secondary vorticities to decrease except where the resulting 
increased path length offset this effect. 

6. The loss coefficient measured in an experiment wth a 
rotating S-shaped pipe showed the same trend, that is, rapid 
increase as the relative velocity ratio W decreased, ns evi- 
denced by the secondary vorticity parameter in the theoreti- 
cal results. 

Le-wis Flight Propulsion Laboratory 

National Ad'visory CoMAnTTSB for Aeronautics 
Cle-ybland, Ohio, August SI, 195S 
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APPENDIX A 


SYMBOLS 

The following symbols are used in this report: | k 

H Bemoulh constant relative to rotating system, 

dimensionless, eq. (4) 
i unit vector in Ui -direction 

unit vector in 'lij-direction 
unit vector in 'ifs-direction 
principal normal of streamline 
total pressure of fluid in stationary channel 
static pressure, dimensional 
radial length, dimensionless 
radial length, dimensional 

, Ui, TJz square roots of components of metric tensor 
associated with orthogonal curvilinear coordi- 
nate sjrstem, eq. (13) 

-Ui, Mj orthogonal curvilinear coordinates, eq. (18) 
fluid velocity in stationary channel 
W velocity relative to rotating system, dimension- 

less 

Z axial distance, dimensionless 

a angle between tangent to streamline and 

direction of axis of rotation 

/3 angle between radial line and streamline, posi- 

tive in direction indicated by positive w 
f total vorticity in stationary channel 

0 coordinate in right-handed cylindrical coordi- 

nate system 

M central angle subtended by streamtube width 


cmvature of streamline 
geodesic curvature of streamline 

X loss coefl&cient 

{ vorticity component parallel to streamline, 

dimensionless 

p fluid mass density, dimensional 

T , ratio of streamtube height in axial direction at 

any point to that at inlet 

<p angle between principal normal of streamline 

and normal to Bemoulh surface,'positive when 
indicating positive rotation about stre amlin e 
according to right-hand rule 
angle between vortex line and streamline in 
stationary channel 

0 total relative vorticity, dimensionless 

CO angular velocity of system, dimensional 

V vector operator 

Subscripts: 

n index equal to 1, 2, 3, denoting cbmponents 

in Ui-, v^, lis-directions, respectively 

1 inlet 

t characteristic point ia channel which is tip in 

case of rotating radial channel 

X ai'biti'ary point along streamtube 

Superscript: 

— vector quantity 
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APPENDIX B 


EXPANSION OF 


V- 


TTXVH 

T7> 


In order to simplify equation (8) and bring out its physical 
significance, the term V • is expanded as follows: 


V • . (R^X VB)+^ VH . fi (Bl) 

This expression is further simplified by expressing the first 
term on the right as 

(1RXVS)=:^[(R7.vW)-6^^XVB)+(PFXQ)-(FFXVS)] 

Since the streamhnes lie on surfaces of constant H (see fig. 2), 
VJ? is normal to W and the preceding equation becomes 

• (Wx VH)=^ [^. . (^X VB) +TR*VB • fl] (B2) 

Equations' (1) and (3) and the fact that VB is normal to W 
are used to simplify the term vB • as follows: 


VB-£2=VB- 


Aw- 




XR^+VBxR^ 


W-W 


Combining equations (Bl), (B2), and (B3) yields 


.-2VB.- 

(I) 

(B3) 


V-^^^P^=^(W"-VW)-a^XVB)+^VB.i (B4) 


The first term on the right of equation (B4) involves a triple 
scalar product in which the term ("PF-vW) represents the 
relative acceleration of the fiuid particle. This acceleration 
has two components, one tangent to the streamline and the 


other along the principal normal equal to W^k, whore k is 
the curvature of the streamline. Hence, the triple product 
is equal to —W^\VH\k sin <p, where <p is the angle measmod 
from the principal normal of the streamline to the normal 
to the BemouUi surface (in the direction of the gradient of 
B), positive when indicating a positive rotation about the 
streamline according to the right-hand rule. The quantity 
K sin ^ is the geodesic curvature k, of the streamline, so that 
this term is zero when the streamlines are geodesics on the 
Bernoulli surface. Hence, equation (Bl) becomes 


WXVB 2 , 2k,\VH\ 


( 9 ) 


Equation (9) is used in the simplification of equation (8) to 
clarify the physical significance of the various terms of 
equation (8). 
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